Population imbalanced Fermi gases in quasi two dimensions 
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We study s-wave pairing of population imbalanced Fermi atoms in quasi two dimensions using 
a mean field theory. At zero temperature, we map out the phase diagram in the entire Bardeen, 
Cooper and Schrieffer-Bose Einstein condensation (BCS-BEC) crossover region by investigating the 
effect of weak atom tunneling between layers. We find that the superfluid phase stabilizes as one 
decreases the atom tunneling between layers. This allows one to control the superfluid-normal first 
order phase transition by tuning a single experimental parameter. Further, we find that a tunneling 
induced polarized superfluid phase appears in a narrow parameter region in the BEC regime. At 
Finite temperatures, we use a Landau-Ginzberg functional approach to investigate the possibility of 
spatially inhomogeneous Fulde-Ferrel-Larkin-Ovchinnikov (FFLO) phase in the weakly interacting 
BCS limit near the tricritical point of spatially homogenous superfluid, FFLO, and normal phases. 
We find that the normal-FFLO phase transition is first order transition as opposed to the continues 
transition predicted in zero temperature theories. 



I. INTRODUCTION 

Since its realization in dilute atomic gases, superfluid- 
ity of alkali atoms has been studied extensively in various 
externally controllable environments For the case of 
fermionic systems, superfluidity arises due to the Bose 
condensation of pairs of fermions at low temperatures. 
One of the most fascinating control parameters has been 
the two-body scattering length between two atoms in two 
hyperfme spin states. The two-body scattering length in 
three dimensions (3D) can be controlled dramatically by 
the use of magnetically tuned Feshbach resonance [2|. 
For dilute atomic systems at low temperatures, the two- 
body interaction is linearly proportional to the scatter- 
ing length in 3D. As a result of this proportionality, by 
controlling the two-body scattering length, the smooth 
crossover between the Bose-Einstein condensation (BEC) 
of strongly bound diatomic molecules to the BCS limit 
of weakly bound cooper pairs has been observed exper- 
imentally. For a negative scattering length, the inter- 
action between two atoms in two hyperfme spin states 
is attractive and momentum space paring gives BCS su- 
perfluidity at low temperatures. For a positive scattering 
length, the interaction is repulsive and two body bound 
states exist in the vacuum gives composite bosonic nature 
for two atoms paired in coordinate space. Bose Einstein 
condensation of these composite bosons gives superfluid- 
ity on the BEC side of the resonance. These two regimes 
smoothly connect at unitarity where the scattering length 
is infinite in 3D. For attractive two body potentials in two 
dimensions (2D), there always exists a two body bound 
state. The 2D bound state energy depends on both 3D 
scattering length and the laser intensity which used to 
create one dimensional lattice to accommodate 2D lay- 
ers. Therefore, 2D paring interactions can be controlled 
by tuning either the 3D scattering length or the laser 
intensity. 

The most recent experiments with ultra-cold Fermi 



gasses have been the focuss of population imbalance 
which leads to the competition between superfluidity and 
magnetism 0, Because of the large spin relaxation 
time of the atoms, experimentalists were able to main- 
tain a fixed polarization P = (iV-j — Ni)/(N^ + 2VjJ over 
the entire time of the experiments, here JVj / j is the num- 
ber of atoms in up/down pseudo spin state. 

For the systems of atoms trapped in external harmonic 
potentials in three-dimensions, phase separation between 
normal phases and various superfluid phases has been 
experimentally observed [H, G|. Theoretical investiga- 
tion of population imbalanced fermion paring in three- 
dimensions (3D) has been encouraged by these series of 
recent experiments For spatially homogenous sys- 
tems, various exotic phases, such as Sarma phase, Fulde- 
Ferrel-Larkin-Ovchinnikov (FFLO) phase, a phase with 
deformed Fermi surfaces, and phase separation have been 
suggested [6J. In trapped systems, various phases are 
separated into concentric shells and the shell structure 
depends on both the interaction strength and the polar- 
ization 0, [H . The boundary between normal and super- 
fluid regions depends on the trap geometry. Experiments 
done in high aspect ratio traps show that the superfluid- 
normal boundary does not follow the equipotential con- 
tours of the trap and show significant distortion of the 
central superfluid region [3j. Quit remarkably, this dis- 
tortion of the central superfluid shell can be explained 
by the surface tension between superfluid and normal 
regions U, EH. Local microscopic physics on 

the superfluid-normal boundary (due to the energy cost) 
causes this surface tension. The most recent theoretical 
studies in 3D reveal the importance of the interaction 
in the normal phase to correctly explain the experimen- 
tally observed Chandrasekhar-Clogston limit of critical 
polarization 13j, |l4| . 

The recent theoretical efforts of understanding the 
fermion pairing in two-dimensional population imbal- 
anced systems attempt to explore the phase diagram in 
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the BCS-BEC crossover region [H, [l7| ■ Experimen- 
tally, such a system can be created by applying a rela- 
tively strong one-dimensional (ID) optical potential to 
an ordinary three-dimensional system. Previous theoret- 
ical studies have concentrated on the nature of various 
phases in 2D layers where the tunneling between layers 
are neglected. 

As the experimental setup in 2D is created out of an 
ordinary 3D system by applying a ID lattice, atom tun- 
neling between layers are always present. This inter-layer 
atomic tunneling can be controlled by a single parame- 
ter, namely the intensity of the ID optical lattice. The 
purpose of the present work is; (1) understand the ef- 
fect of tunneling, (2) find out how one can control the 
phase transitions by using this easy controllability, and 
(3) study possible FFLO phase in 2D. The motivation 
of studying FFLO phase came from the fact that FFLO 
phase is more favorable in low dimensions (l8| . In 3D, 
FFLO phase has not been seen in experiments, though 
theory predicts that such a phase can exist in a narrow 
parameter region in the weakly interacting limit. 

The first part of this work is a natural generalization of 
the work presented in refs. [H, [3 E3] t° include the weak 
atom tunneling between atomic layers. We find that the 
first order superfluid-normal phase transition can be ef- 
fectively controlled by the laser intensity. In typical 3D 
population imbalanced gases, one has to change the pop- 
ulation imbalance to control the superfluid-normal phase 
transition. Moreover, we find that the inclusion of tun- 
neling allows to stabilize the polarized superfiuid phase 
in BEC regime. Further, in the weak coupling BCS limit, 
we investigate the possible inhomogeneous FFLO phase 
near the tricritical point of superfiuid, FFLO, and normal 
phases and find that the normal-FFLO phase transition 
is of first order. 

The paper is organized as follows. In section II, we 
consider a zero temperature mean field theory and de- 
rive analytical expressions for the energy, density and 
gap equation for weakly coupled layers. Then we predict 
the phase diagram in the entire BCS-BEC crossover re- 
gion and discuss the effect of weak atom tunneling on the 
superfluid-normal boundary. In section III, we neglect 
the coupling between layers and use a Landau- Ginzberg 
functional approach at finite temperatures. In the weak 
coupling BCS limit, we derive an analytical expression 
for the Landau's free energy functional and discuss the 
possible FFLO phase near the tricritical point. Finally, 
our conclusion is given in section IV. 

II. ZERO TEMPERATURE PHASE DIAGRAM 

We consider an interacting two-component Fermi 
atomic gas trapped in quasi two dimensions. The ID op- 
tical potential created by the counter propagating laser 
beams has the form V — sEnsiv?(2'Kz/X). Here A is the 



wavelength of the laser beam, Er = h 2 (2ir / A) 2 /(2M) 
is the recoil energy. The dimensionless parameter s can 
be used to modulate the laser intensity. When the pa- 
rameter s is large, the atomic system forms a stack of 
weakly coupled 2D planes with periodicity d = A/2. The 
Hamiltonian of the system H = Y^j Hj is represented by, 

H i = / dV {£4^)h^f -AOMr) (1) 

^ (7 

+*X^4,( r )V'i+ia(r)+/ic] 

+t/ 2 £>V'jt M^jl ( r )^'i M^'T ( r )| 

where j is the layer index with r 2 = x 2 +y 2 , V2D is the 2D 
gradient operator and U2D is the 2D interaction strength. 
The operator (r) creates a fermion of mass M in j'th 
plane with pseudo spin a =T,j at position r = (x,y). 
We consider a tight ID lattice where the atomic wave 
function becomes more and more localized in the planes. 
Using the harmonic approximation around the minima of 
the optical lattice potential [13], we find the lattice tun- 
neling parameter t/En — (2s 3 / 4 / y/n) exp[— 2y/s\. Notice 
that the inter layer tunneling energy t can be varied by 
changing the laser intensity parameter s. In the limit 
t — ► 0, the system is decoupled planes of Fermi atoms. 

In this section, we consider zero temperature and use a 
mean filed theory to decouple the interaction term writ- 
ing U2D^}- l (r)tl>} l {r)^ j i(r)i/} j f(r) = [A^ n (r)^(r) + 
h.c] — \A\ 2 /U~2D- Further we neglect the possible inho- 
mogeneous Fulde-Ferrell-Larkin-Ovchinnikov paring and 
leave the finite temperature FFLO discussion for the next 
section. The Fourier transform of the Hamiltonian gives, 

H= ( Cfe ll _ ^) a lna( k \\) a "^{k\\) (2) 

feii ,er,m 

(k\\)a mrT (k\\) + h.c] 

fcy ,£T,m 

+ J2[Aal t (k ll )al l (-k ll )+h.c]- 1 -^- 

tey ,m 

where e fe|| = ft 2 fc 2 /(2M) with fc| = fc 2 +fc 2 . The notations 
m and fey represent the index of the plane and the mo- 
mentum parallel to the 2D layers. The operator a] ri(7 (fe||) 
creates a fermion of mass M in mth plane with pseudo 
spin a =Tii an d momentum fcii. The periodicity along 
the z-direction allows us to write the Fermi operators, 

a mT(fc||) — exp[ik z md]c^(k) 

a ml(-fc||) = ^exp[-ifc z md]q(-fc) (3) 

fc* 
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FIG. 1: Phase diagram of quasi 2D Fermi system in (a) BCS 
regime, (b) Unitarity regime, and (c) BEC regime, where we 
fixed the 3D scattering length to be a s = — \Jti 2 /{2mE R ), 

a s — ► oo, and a s = + 'ft 2 / '(2mE 'r) respectively. The lattice 
height is fixed by setting s — 8. The abbreviations are; SF: 
superfluid phase, PSF: polarized superfluid phase, NM: nor- 
mal mixed phase, PN: fully polarized normal phase, and O: 
vacuum state. 



where the fermi operator cj.(— k) creates a fermion of 
momentum k = (k x ,k y ,k z ) with pseudo spin a. Trans- 
forming the Hamiltonian given in Eq. ^ using above 
canonical transformation followed by the usual Bogoli- 
ubov transformation, the Hamiltonian per plane can be 
expressed in terms of quasi particle Fermi operators at CT 
as 



H/N= E R 
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where N is the number of planes, e k 
and E k ± - 



e k + 2tcos(k z d), 



70 
60 
« 50 

EH 

* 40 
30 
20 



(b) 



SF 



NM 



FIG. 2: Chemical potential on the superfluid phase and nor- 
mal phase boundaries as a function of lattice height s. We 
fixed the parameters in the BCS regime a s = 
Panel (a) shows the chemical potential on the boundary be- 
tween superfluid phase and the fully polarized normal phase 
at h — Er and panel (b) shows the chemical potential on the 
boundary between superfluid phase and the normal mixed 
phase at h = AEr. We find same qualitative behavior in the 
entire BCS-BEC regime. 



potential fi = (/x-j- + Ml)/2 and chemical potential differ- 
ence h = (/if — /ij.)/2. Without loss of generality, we take 
h > 0. The energy density per plane — (— l/flo) M-^g] 
with Z G = Trexp[-f3 H/N] is 



2^ J (2^ 2 

-T 



h - M - V(efe - M) 2 + A 2 



(5) 



c/ 2D A) 



d 2 k l{ 



ln[l + exp(-/3 S fe+ )] 
+ ln[l + cxp(/3 ^fe-)] 



Here /3o = ^-/{ksT) is the inverse temperature. Due to 
the nature of the contact interaction we used, Eq. ([S]) 
is divergent so that proper regularization must be done. 
This can be done by relating the 2D contact interaction 
U2D to the bound state energy Eb as [2(| 



U 2D 



d 2 k 



(2tt) 2 % 2 k 2 /M + Eb 



(6) 



/i ± a/ (efc — /i) 2 + A 2 . Average chemical 



The bound state energy is given by Eb = 
{CTiijJl/'k) exp[\/27r'L/ a s], where a s is the 3D s- 
wave scattering length, uj^ — \J^ 2 sEr] (m\ 2 ) is the 
trapping frequency due to the lattice, II = Ti / {muj l) 
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is the oscillator length and C w 0.915 211 ) . Notice that 
the bound state energy depends not only on the 3D 
scattering length but also on the lattice potential. After 
replacing 2D contact interaction Uid in Eq. ([5]) by Eq. 
and then combining the first and second term in Eq. |(SJ), 
we remove the ultraviolet divergency. Converting the 
sum over k z into an integral over the first brillouin zone 
and then performing the both momentum integrals, the 
energy density at zero temperature is 



2nh z 



2 4 2 V ^ 



(7) 



Q(h - A)hVh 2 - A 2 
3AV + 2/T 3 



2(M 2 
15AV 

"8(^ 2 + A 2 ) 7 /2 



A 2 ) 3 / 2 
4 + G(i 6 



where Heaviside theta function Q(x) — 1 for x > and 
otherwise. Then the gap equation, the number density 
and density difference are calculated by using the equa- 
tions, <9fi/9A = 0, n = —dfl/dfi, and = —dCl/dh 
respectively. 
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(9) 



^e(^A)\//i 2 -A 2 



(10) 



We solve Eq. (|8|) for given s, a s , fi and ft for A and 
then determine the energy of the system from Eq. 0. 
Comparing the energy of the system, we then map out 
the phase diagram in /i — h parameter space as shown in 
FIG. Q] As can be seen in FIG. HJ the qualitative feature 
of the phase diagrams in the entire BCS-BEC regime is 
the same. However the normal mixed phase (where both 
spin up and down components coexist) narrows down as 
one goes from BCS regime to BEC regime. In trapped 
systems, h is uniform within the cloud, while chemical 
potential \i varies monotonically from the center to edge 
of the trap. Taking vertical slices of the phase diagram at 
fixed h, one can see that the atomic cloud passes through 
various local phases giving a shell structure in a trapped 




FIG. 3: The density of a quasi 2D Fermi system as a 
function of chemical potential fi in (a) BCS regime [a* = 
-•\A 2 /(2mE R ), and h = 2E R ] and (b) BEC regime [a* = 
+ h 2 / (2mE r) , and h = 6Er]. The lattice height is fixed by 
setting s = 8. 



system. The shell structure depends on both interaction 
strength Eb and the chemical potential difference h. The 
similar phase diagram in the fi—h plane is obtained in the 
limit of s — > oo in Ref . [l5( and Ref. [l6( . The main differ- 
ence is the superfluid phase region in the phase diagram 
narrows down as a result of weak tunneling (due to the 
finite s) . In contrast to 3D population imbalanced Fermi 
systems 0], the parameter region for polarized super- 
fluid phase (both superfluid and normal coexisted phase) 
in the BEC regime is very small. Using a very similar 
mean field analysis, Ref. [l5( and Ref. [l6( conclude that 
the polarized superfluid phase is absent in the limit of 
s — > oo in 2D. However, using a beyond mean field ap- 
proach, Ref. 17| predicts a polarized superfluid phase 
oo limit. Nevertheless, within our mean 



Ref. 

even in the s 

field description, we find a tunneling induced polarized 
superfluid phase only in the BEC regime. Notice that 
we have constructed the phase diagram by solving the 
equations for fixed /i and h. Alternatively, one can use 
the canonical ensemble (where n and are held fixed) 
to construct the phase diagram in n-rid plane. As is 
zero in the SF phase, a line in n-rid plane will represent 
the SF phase. Between the SF line and PN/NM phase 
separates a region of phase separation [16( . 

Figure [5] shows the variation of the chemical poten- 
tial on the superfluid-normal boundaries. As the lattice 
height increases, the superfluid phase stabilizes against 
the normal phases. As a result, the superfluid cloud ex- 
tends toward the edge as one increases the ID lattice 
height in a trapped system. This is because, 2D inter- 
action strength increases with increasing the ID lattice 
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FIG. 4: The coefficients a (solid line), 7 (dashed line), and v 
(dotted line) as a function of temperature. We fixed the value 
of chemical potential difference h — 1.08fcsTb such that there 
are three phase transitions as lower the temperature. See also 
FIG. E] 

height. As the lattice height is controllable though the 
laser intensity, the laser intensity can be considered as 
a non-destructive experimental knob to control the first 
order superfluid-normal phase transitions as well as po- 
larized supcrfluid phase. This easy controllability is avail- 
able only in quasi 2D systems and one has to change the 
population imbalance in 3D systems to control the tran- 
sition. 

In figure [31 we present typical density profiles for two 
different representative values of interaction strengths in 
the BCS regime and BEC regime. For qualitative un- 
derstanding, the chemical potential axis can be consid- 
ered as a spatial coordinate in trapped systems. This is 
because, as we mentioned before, the chemical potential 
varies monotonically from the center to edge of the cloud. 
The density profile in the BCS regime shows a kink and 
a discontinuity representing two phase boundaries. The 
discontinuity represents the phase boundary between su- 
pcrfluid phase and the mixed normal phase while the kink 
represents the phase boundary between mixed normal 
phase and the fully polarized normal phase. In contrast, 
the density profile in the BEC regime shows a disconti- 
nuity showing a two-shell structure at given h and s. 

III. FFLO PHASE NEAR TRICRITICAL POINT 

In this section, we consider a finite temperature spin 
imbalanced Fermi system in two dimensions and neglect 
atom tunneling between layers. The Hamiltonian of the 
system in the mean field description reads, 

H = j ^{E^Wt-Hf -lAMr) (11) 

+<rhil)l{r)%l)< T (r) + A{r)i/j\(r)ip\(r) + ft.c| 
where a is the pseudo spin f (+) and \ (— ). In order to 



study the FFLO state near the tricritical point, we use 
a Landau's phenomenological approach to write the free 
energy functional in the weakly interacting BCS limit. As 
shown in Ref. [2fJ], the chemical potential of a 2D Fermi 
gas is given by [i = ep — Eb/2, where tp and Eb are 
Fermi energy and binding energy of the two body bound 
state. In the weak coupling BCS limit, Eb <C ep, so that 
the chemical potential can be approximated by the Fermi 
energy. Further, spatial modulation of the order parame- 
ter A(r) is small close to the tricritical point. Therefore, 
free energy can be expanded in powers of the superfluid 
order parameter. The coefficients of the terms in each 
order determine the nature of the each phase transition. 
Following the reference (2^ |. the thermodynamic poten- 
tial can be written as, 

F = a|A| 2 + 7 ^|dA| 2 /2 + 7 |A| 4 ( 12 ) 
+u{\A\ 6 + 3vp\d 2 A\ 2 /l6 + 2^|A| 2 |aA| 2 
+v 2 p[(A^ 2 (dA) 2 + A 2 (dAt) 2 ]/4} 

with the coefficients a = N(0){\n(T/T c ) + Reip[l/2 + 
ih/(2itT)] - ^(1/2)}, 7 = N(0)irK 3 /4 and v = 
— N(Q)ttK^/8. Here ip is the digamma function, T is 
the temperature, T c is the critical temperature, 1V(0) is 
the density of states at the Fermi surface, vf is the Fermi 
velocity, and the functions K3 and K5 are given by, 

Kz = -^^Re[^ 2 \x)] (13) 

K 5 = ±—.R e [ 1 hW(x)] 

Here x = 1/2 - ih/(2Trk B T) and rp^(x) is the n-th 
derivative of digamma function on its argument. No- 
tice that the forth order term 7 simultaneously accom- 
panied by the gradient term (|<9A| 2 ) in the BCS mean 
field theory so that we must include the sixth order term 
v. When both a and 7 are positive, thermodynamic po- 
tential gives a single minimum at A = and the sys- 
tem is in the normal state. We consider the simplest 
FFLO state where the Cooper pairs in the coordinate 
space have the plane wave form A(r) = Aoexp[ig • r\. 
When 7 becomes negative, the modulated order has the 
lowest energy. The tricritical point is determined by the 
conditions a = and 7 = 0. The second order phase 
transition from normal state to homogenous superfluid 
state is determined by the condition a = 0. Above 
the superfluid transition (T > Tc), a > and below 
the transition a < (see FIG. UJand discussion below). 
For the first order homogenous superfluid-normal phase 
transition, we find the value of A at the transition by 
setting dF/dA — at q — 0. This leads to the solu- 
tion A± = [— 7 ± \/ j 2 — 3av]/(3v). The condition for 
the first order transition is then determined by setting 
F(A = 0) = F(A + ) which gives 7 = -2^Jav. 
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FIG. 5: Phase diagram of 2D population imbalance Fermi 
system near the tricritical point. All three phases, superfluid 
(SF), normal, and inhomogeneous FFLO coexist at the tri- 
critical point represented by the solid circle. 



Let us now consider the inhomogeneous superfluid- 
normal phase transition with finite q. Using above simple 
anzats for the order parameter A, the free energy is given 
by F = (a + 7 Q 2 /2 + 3z/Q 4 /16)A 2 + (7 + 3^Q 2 /2)A 4 + 
1/A 6 , where Q = vfQ- For a possible second order phase 
transition, we set the coefficient of A 2 to be zero, and 
then minimize it with respect to Q. This gives the condi- 
tion for possible second order transition into FFLO state; 
7 = — 2-\/3<w/4 with center of mass paring momentum 

q = v /-8 7 /(6i/)/'Uf- 

For a possible first order transition into FFLO state, 
we minimize the free energy with respect to both A and 
Q. This gives the condition 7 = —2^/6^/5 with paring 
momentum q = \f— 2j/(6v)/vf- 

Transforming these second order and first order phase 
transition conditions into temperature (T) and chemical 
potential difference (h) through T — h dependence of the 
parameters, a, 7, and v, we find that the first order phase 
transition line is always appears before the second order 
phase transition line in T — h plane as one decreases h. 
This indicates that the transition from normal phase into 
FFLO phase in a 2D population imbalanced system at fi- 
nite temperature is first order. This contrast to the zero 
temperature theories [lCil l2~ij , where thermodynamic po- 
tential is expanded only to the second order in A. In 
these zero temperature theories, the transition from nor- 
mal state into FFLO state is found to be a second order 
in 2D. Notice, here we used the simplest FFLO state 
which has the plane wave form. We do not expect that 
the nature of the transition would change if we choose 
somewhat complex forms of the FFLO state. 

The phase diagram in the vicinity of tricritical point 
is shown in FIG. [5j As our theory is not valid at low 
temperatures, we present our results at a finite temper- 
ature range. In FIG.JH we show the temperature depen- 
dence of the parameters a, 7, and v. As these parameters 
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FIG. 6: Center of mass paring momentum as a function of 
temperature for h/(kBTc) = 1.11 (dotted line), /i/(fesTcr) = 
1.10 (dashed line), and h/(k B T c ) = 1.09 (solid line). 

are functions of both temperature and chemical potential 
difference h, we choose h(= l.08kgTc) such that there 
are three phase transitions as one lowers the tempera- 
ture. The sequence of these phase transitions are Normal 
— > FFLO — > SF — > FFLO. The transition temperatures 
for Normal-FFLO, FFLO-SF and SF-FFLO are 0.538T C , 
0.533Tc, and 0. 434Tb respectively. Notice that the pa- 
rameter a is negative in the SF phase and positive in 
both normal and FFLO phases. Further, 7 changes its 
sign as the system undergos first order transition into 
FFLO state. The center of mass pairing momentum in 
the FFLO phase is shown in FIG.[6]for three different rep- 
resentative values of chemical potential differences. As 
shown in the phase diagram in FIG. [51 the FFLO window 
always appears in between superfluid and normal phases. 
Therefore, in trapped systems, the FFLO phase will ap- 
pear bordering between superfluid and normal concentric 
shells at low temperatures below tricritical point. 



IV. CONCLUSIONS 

We studied Fermi superfluidity in quasi two dimen- 
sions focussing on population imbalance and weak atom 
tunneling between layers. By using a mean field theory, 
we derived an analytical expressions for the free energy, 
density, and solution of the gap equations at zero tem- 
perature. Analyzing the free energy of the system, we 
mapped out the phase diagram where we find phase sep- 
aration due to the population imbalance. Further, we 
find that the suppression of weak atomic tunneling sta- 
bilizes the superfluid phase due to the enhancement of 
the 2D interaction strength. The easy controllability of 
the tunneling through the laser intensity allows one to 
tune the first order supcrfluid-normal phase transition. 
We also find that the polarized superfluid phase which is 
absent in pure 2D limit in the mean field description, can 
be realized by allowing atoms to tunnel between layers 
in the BEC limit. The tunneling dynamic of the system 
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can be understood by generalizing the results in refer- 
ence 

At finite temperature, we used a Landau's functional 
approach in the weak coupling BCS limit and dis- 
cussed the inhomogeneous FFLO phase near the tricrit- 
ical point. We find that as one decreases the population 
imbalance, system undergoes a first order phase transi- 
tion into FFLO phase where zero temperature theories 
predict a second order phase transition. 
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